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Abstract 

Results of an experimental search for new macroscopic forces with Yukawa range between 5 and 
500 microns are presented. The experiment uses 1 kHz mechanical oscillators as test masses with 
a stiff conducting shield between them to suppress backgrounds. No signal is observed above the 
instrumental thermal noise after 22 hours of integration time. These results provide the strongest 
limits to date between 10 and 100 microns, improve on previous limits by as much as three orders 
of magnitude, and rule out half of the remaining parameter space for predictions of string-inspired 
models with low-energy supersymmetry breaking. New forces of four times gravitational strength 
or greater are excluded at the 95% confidence level for interaction ranges between 200 and 500 
microns. 
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I. INTRODUCTION 



Experimental tests of Newtonian gravity and searches for new weak forces in addition to 
gravity have been conducted over length scales ranging from light-years down to laboratory 
distances. Based on recent results fl] and reviews [2J], new forces with strength weaker than 
or comparable to gravity have been excluded over distances ranging between 200 /im and a 
light-year. To date, only the single experiment in Ref. [l| has attained gravitational sensitiv- 
ity below 1 mm, and limits on the strength of new interactions increase very rapidly below 
100 /im |3(. The sub-millimeter range attracts continuing experimental interest because 
many recent theoretical attempts at the unification of fundamental forces predict specific 
new phenomena in this regime. We present the results of an experiment which has attained 
a maximum sensitivity of about four times gravitational strength above 200 /im, and which 
provides the most sensitive limits on new forces between 10 and 100 microns. 

A. Current Limits 

Results from experimental searches for new macroscopic forces are most commonly param- 
eterized by a Yukawa interaction. The potential due to gravity and an additional Yukawa- 
type force between two mass densities pi(rl) and ^2(^*2) can be written 

Gpi(f 1 )p 2 (f 2 ) 

where G is the Gravitational constant, 7*12 is the distance between f± and r*2, a is the strength 
of the Yukawa interaction relative to gravity, and A is the range. Current experimental limits 
on a between 1 /im and 1 cm are shown in Fig. [7J together with several recent theoretical 
predictions of new phenomena. 

All previously published limits in Fig. UJ (bold, dashed curves) are obtained from torsion 
balance experiments. These curves represent 95% confidence level limits with the exception 
of the result from the Irvine 2-5 cm null experiment j^], which is a 1 a limit. The torsion 
balance experiment at the University of Washington has attained gravitational sensitivity 
or better for ranges above 200 /im [l|. Below 20 /im the best previous limit is derived from 
the Casimir force measurement by Lamoreaux 51; the curve shown in Fig. U\ is based on the 

no 

analysis by the authors p, |]| • From the figure, the previous experimental limits allow for 
forces in nature several million times stronger than gravity at ranges as large as 20 /im. 
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B. Motivation 



In addition to the unexplored parameter space, there are many specific motivations to 
search for new effects below 100 /im Several arise from the decades-long effort to 
describe gravity and the other fundamental interactions in a unified theoretical framework. 
The leading candidates for such a theory are string or M theories, which must be formulated 
in more than three spatial dimensions. In the traditional view, the extra dimensions are 
taken to be compactified on a scale on the order of the Planck length (Lp ~ 10 -35 m 
or Lp 1 « 10 19 GeV), a fundamental scale at which the unification is expected to occur. 
Recently, a class of models has been discovered in which the unification can occur at much 
lower energies, near the weak scale (L^ 1 « 1 TeV) jjj. In this picture, perhaps the most 
frequently cited development in particle physics in recent years, the discrepancy between 
the Planck scale and the weak scale (the so-called hierarchy problem) is removed as a 
consequence of some of the extra dimensions remaining large and accessible only to gravity, 
with the Standard Model fields confined to the usual three dimensions Q]. The size R 
of the extra dimensions is given by R n = L^y n /Lp, where n is the number of compact 
dimensions. For n = 1, R ~ 10 13 m, clearly ruled out by astrophysics. However, the choice 
of n = 2 implies R ~ 1 mm, with the consequence that the gravitational potential will 
behave according to a 1/r 3 law below this scale. As the scale R is approached from above, 
Yukawa corrections are predicted; the model illustrated in Fig. [7| predicts a = 4 [ill ]. 

Another class of predictions arises in superstring theories in which supersymmetry 
(SUSY) is broken at low energies. These models are of interest as unification scenarios 
that account for the absence of flavor-changing neutral currents. Superstring theories gen- 
erally contain gravitationally-coupled scalar fields called moduli, which are massless at the 
string scale but which can acquire mass from the same process which breaks SUSY. For 
models in which SUSY is broken by a gauge-mediated process between 10 and 100 TeV, 
the moduli acquire masses at the sub-eV level corresponding to interaction ranges between 
100 /im and 10 cm, with values of a as high as 10 6 [12]. Predictions for the largest such 
effects (for the moduli which couple to the strange quark and gluon) with the uncertainties 
are shown in Fig. [7| In a related model (radius modulus) in which SUSY is broken near 
1 TeV via weak-scale compactification 13j] , the moduli can acquire Compton wavelengths 
in the range from 10 jum to 1 mm. The associated strength is predicted to be a = 1/3, as 
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shown in Fig. [7| 

The dilaton is another scalar field predicted by string theories, which can also acquire 
mass from SUSY breaking. Due to its more universal couplings, it is expected to acquire 
a mass too large to be observable in the scenario in Ref. 12j, but other authors take the 
dilaton mass to be an unknown parameter. The estimated coupling for the dilaton from 
Ref. [14] is shown in Fig. [7J where previous experimental limits exclude a mass less than 
about 3.2 x 10~ 3 eV; more recent calculations can be found in Ref. jis| . 

The axion, a light pseudoscalar boson motivated by the strong CP problem of the Stan- 
dard Model, is also predicted to mediate macroscopic forces. Astrophysical and laboratory 
bounds have left an allowed window for the axion mass corresponding to 200 /im < A < 
20 cm Q|. The coupling of the associated long-range force between unpolarized test masses 
is constrained by measurements of the neutron electric dipole moment to the region indicated 
in Fig. For pseudoscalars which are not derivative coupled (unlike the axion), double ex- 
change leads to a 1/r 3 potential between unpolarized test masses; short range experiments 
may set better limits on pseudoscalars via this effect |17| . 

Additional predictions are motivated by the cosmological constant (A) problem, the dis- 
crepancy between the observed flatness of the universe versus the extreme curvature expected 
from the vacuum energy contributions of the Standard Model fields. Recent models jisl Q] 
assert that if a sufficient fraction of the energy density in the universe is in the form of 
vacuum energy, consistency with local field theory implies the existence of new interacting 
quanta with mass on the order of A 1//4 . The predictions in Fig.Elare derived from the conser- 
vative assumption of cosmological vacuum energy density p\ in the range 0.1p c < pa < p c , 
where p c is the critical density of the universe. 



II. EXPERIMENT 



A. Apparatus 

The present experiment is illustrated in Fig. ^ The test masses consist of a planar 
resonant detector mass and a planar source mass driven at the detector resonant frequency. 
This geometry is chosen to concentrate as much test mass density as possible at the length 
scale of interest. It is also largely null with respect to 1/r 2 forces and so is effective in 
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Tuning Block 



FIG. 1: Central components of the apparatus 

suppressing the Newtonian background in the context of a new force search. The source 
mass consists of a 35 mm x 7 mm x 0.305 mm node-mounted tungsten reed. It is driven 
at its second cantilever mode by a PZT bimorph at a frequency carefully matched to the 
1 kHz resonant frequency of a normal mode of the detector mass. The detector mass consists 
of a 0.195 mm thick tungsten torsional oscillator of a double-rectangle design. The small 
rectangle measures 11.455 mm x 5.080 mm. For the resonant mode of interest (the 5th 
normal mode) the two rectangular sections counter-rotate about the torsional axis, with 
most of the amplitude confined to the motion of the small rectangle. This double-torsional 
arrangement provides si gnifi cant isolation of the small rectangle from the detector mount, 
reducing mode damping |20j|. When data are collected, the front edge of the source mass is 
aligned with the back edge of the small detector rectangle, and a long edge of the source is 
aligned above the detector torsion axis. This geometry (called "fiducial" below) maximizes 
the on-resonance torque on the detector if a mass-coupled force is present between the 
source and detector. 

A stiff conducting shield suspended between the test masses suppresses electrostatic and 
acoustic backgrounds. The shield consists of a 0.060 mm thick sapphire plate with a 100 nm 
gold plating. The shield is mounted on a brass support and secured on opposite ends with 
epoxy. 
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The test masses and shield support are each suspended from modified optical mirror 
mounts, which are used as a tilt stages to level each element in the horizontal plane. Each 
tilt stage is in turn mounted to the bottom segment of a vibration isolation stack. Each 
stack consists of five solid brass disks connected by fine steel wires, as described in detail 
elsewhere j^]. At the 1 kHz operational frequency of the experiment, each stack provides 
an attenuation of approximately 200 dB. The vibration isolation stacks of the source and 
detector masses are suspended from inverted 3-axis micrometer stages, which provide full 
translation control. The relative positions and coplanarity of the test masses and shield 
are ascertained by touching the elements against each other, or against a series of 0.5 mm 
diameter sapphire hemispheres attached to small rods at the end of the shield support. 

To further reduce acoustic backgrounds, the apparatus is placed in a 75 liter vacuum bell 
jar, which is pumped to a pressure of approximately 2 x 10~ 7 torr with a liquid nitrogen- 
trapped diffusion pump. In order to isolate the (non-magnetic) central apparatus from stray 
fields generated in the steel pump components, the bell jar is connected to the diffusion pump 
via aim long, 10 cm diameter aluminum riser. The position-control micrometer stages are 
manually operated with torque rods which exit the bell jar via rotary feedthroughs. 

In the absence of electromagnetic, acoustic, and vibrational backgrounds, the experiment 
is limited by thermal noise due to dissipation in the detector mass. To reduce this dissipation, 
the detector mass is annealed at 1300 °C for several hours in an induction furnace under 
a helium atmosphere before installation in the bell jar. This is observed to increase the 
detector mechanical quality factor (Q) typically by factors of 5 or more. 

The detector temperature is stabilized with an electronic temperature controller using 
a silicon diode sensor and resistive heating element mounted to the bottom stage of the 
detector vibration isolation stack. The temperature is maintained at 305 K with fluctuations 
of about 0.1 K. 

B. Readout 

The detector readout electronics and source mass PZT drive are illustrated in Fig. |21 
Oscillations of the detector mass are read out with a capacitive transducer. The capacitive 
transducer probe consists of a 2.5 mm diameter brass cylinder supported with its flat end 
approximately 0.1 mm above a rear corner of the large rectangle of the detector mass. The 
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FIG. 2: Major components of the detector readout and PZT drive 

probe is biased at 200 V through a 100 GQ resistor, whose value must be large to reduce 
current noise. The front end of the preamplifier consists of an SK152 JFET located in a small 
box immediately above the detector oscillator. The 100 mK noise temperature of the SK152 
preamplifier ensures that this circuit is more than sufficient for detecting the 300 K thermal 
oscillations (amplitude ps 100 fm) of the detector mass. The JFET preamplifier is followed 
by a second preamplifier (Stanford Research SR560), filters, and finally a two phase lock-in 
amplifier. The total voltage gain from the capacitive probe to the lock-in input is about 
1600. A crystal-controlled oscillator provides a reference signal for the lock-in amplifier and 
drives the source mass PZT through a 1:10 step-up transformer. 



III. DATA SAMPLE 



A. Force Measurement Data 



The data were recorded in 108 "cycles" over five days. Each cycle consists of 7 runs 
containing both force measurement data and diagnostic data. The diagnostic data are used 
to continually monitor the performance of the instrument. For all runs the bandwidth of the 
lock-in amplifier was set to 250 mHz (roughly 5 times the width of the detector resonance) 
to include the noise power of the detector thermal oscillations, which is used for calibration. 
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A data sampling rate of 1 Hz was used, corresponding to an over-sampling of about a factor 
of 10. 

Within each cycle, the runs were taken in the following order: First, a stable dc bias (usu- 
ally 5-10 V) was applied to the shield to induce a large resonant electrostatic signal. A run 
of 120 samples was taken with the PZT drive frequency set 30 mHz below the approximate 
value of the detector resonant frequency. Four subsequent runs were taken, each with the 
PZT drive frequency incremented by 15 mHz relative to the previous run in order to display 
the detector resonance. The precise value of the detector resonant frequency was determined 
from the five biased runs, and was then used as the drive frequency for the next run of the 
cycle. The shield was then grounded and a run of 720 samples (12 minutes) was recorded. 
This was followed by a shorter run of 288 samples (4.8 minutes) taken with the PZT drive 
set to 1171.000 Hz (at least 2 Hz below the detector resonance), in order to monitor the 
system for non-resonant systematic offsets. The entire cycle was repeated indefinitely until 
an intervention to service the experiment was necessary. A total of 108 such cycles were 
acquired yielding a total of 77760 on-resonance samples. 

A plot of the data from the biased diagnostic runs is shown in Fig. El (only data from 
10 consecutive cycles are shown for clarity). The signals from the two phase-quadrature 
channels of the lock-in amplifier are plotted against each other to show the phase behavior 
of the signal with drive frequency. In all cycles, the signal maximizes at very nearly the same 
phase and magnitude, indicating good stability of the detector resonant frequency, source 
mass amplitude, and system gain. 

Fig. 0] shows histograms of the on- and off-resonance unbiased data. The plots combine 
data from all 108 on-resonance and off-resonance runs, each plot displaying data from a 
single channel of the lock-in amplifier. The data exhibit smooth gaussian behavior centered 
about common means, as expected in the exclusive presence of detector thermal noise and 
amplifier noise. The widths of the distributions of the on-resonance data (left-hand plots 
in Fig. EJ) are roughly twice those of the off-resonance data (right-hand plots) due to the 
contribution of detector thermal noise. The means of the distributions in Fig. H] are shown 
in Fig. in which the data from the separate lock-in amplifier channels are plotted against 
each other as in Fig. El The on-and off-resonance means agree within the 1 a standard 
deviations shown, indicating the absence of any resonant force signal. The offset from the 
origin of the mean of the off-resonance distribution is due to leakage of the reference signal 
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FIG. 3: Data from the biased runs of 10 consecutive cycles of the data sample. Points from 
separate cycles are labeled with unique symbols. 1 a error bars are shown for one cycle. Groups 
of points corresponding to a particular drive frequency are circled. 

internal to the lock-in. 

B. Consistency Checks 

After the acquisition of the session data, several checks are made to verify the stability of 
the results and their consistency with known backgrounds. First, several additional cycles 
are acquired with variations in the test mass geometry, including larger vertical gaps and 
different overlap configurations. No resonant signal is observed in any of these sessions. 
This strongly disfavors the possibility that the observed null result is due to a fortuitous 
cancellation of surface potential, magnetic, and/or acoustic effects, all of which are expected 
to have different dependencies on the geometry. 

Next, several on-resonance runs are acquired with different transducer probe bias voltage 
settings and with the source mass drive turned off. The observed linear dependence of the 
rms fluctuation in these data on the probe bias voltage is consistent with detector motion 
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FIG. 4: Distributions of data samples without shield bias. Left-hand plots: Drive (and lock- 
in reference) tuned to detector resonance. Right-hand plots: Drive tuned 2 Hz below detector 
resonance. 

due only to thermal noise and rules out additional motion from transducer back-action noise. 
This check is important because the magnitude of the detector thermal motion is used for 
calibration. 
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FIG. 5: Means of the distributions in Fig. 0] with channels plotted against each other to show 
phase. Error bars are 1 a standard deviations of the means. The point with the larger error bars 
is the on-resonance mean. 

The data from the diagnostic runs with shield bias voltage applied can be used to estimate 
the minimum size of the residual potential difference between the shield and the (grounded) 
test masses needed to produce a resonant signal. From Fig. |3J the maximum signal with 
9.6 V shield bias is about 1.3 mV. Further measurements show that this signal scales as the 
fourth power of the shield bias voltage, as expected for an electrostatic force between source 
and detector mediated by a deflected shield. Scaling the ratio of the on-resonance error bar 
in Fig. |5]to the 1.3 mV diagnostic implies that a residual shield bias of at least 1.5 V would 
be needed to generate a systematic effect above detector thermal noise. This is about an 
order of magnitude larger than the measured residual potential difference between the shield 
and test masses. 

Magnetic effects can generate background signals though several mechanisms. The most 
important effect involves generation of eddy currents when the source mass moves in an 
external magnetic field. Fields produced by the source eddy currents create eddy currents 
in the detector, which then interact with the applied field. Studies of this effect with 
large applied fields show that the induced force varies as the square of the applied field (as 
expected), and extrapolation to the ambient field actually present indicates that this effect 
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should be about 5 times smaller than the thermal noise limited sensitivity. 
C. Auxiliary Measurements 

Before acquiring the session data samples, several measurements are made to support 
the analysis of the data for evidence of a net force. These include the relative phases of the 
motion of the test masses, the detector mechanical quality factor, and a precise survey of 
the test mass geometry. 

The Yukawa interaction for the test mass geometry of the experiment is used to model 
any potential signal. Such an interaction can result in either a purely attractive or repulsive 
force between the source and detector masses. In order to ascertain the phase of the signal for 
such a force, the source and detector are placed in their fiducial positions with an additional 
vertical offset so that both are situated above the electrostatic shield. An electrostatic bias 
of 1.5 V is applied to the detector mass with a dry cell in order to induce a large, stable, 
purely attractive electrostatic coupling between the test masses. A set of runs is taken with 
small increments of the PZT drive frequency between successive runs. The position of the 
data point in the complex plane corresponding to the maximum signal as a function of drive 
frequency is recorded. The phase of this signal corresponds to that of a purely attractive 
force. The phase of a purely attractive signal is 189°, just below the negative horizontal axis 
in the plots in Figs. El and El 

The detector mechanical quality factor Q is measured by applying a 100 mV resonant 
ac signal directly to the detector, switching the detector to ground, and observing the ring- 
down signal on the lock-in amplifier with the reference frequency tuned 1 Hz off the detector 
resonance. A value for Q is obtained from a least-squares fit to the ring-down waveform. 
Several such measurements are made before and after the data cycles reported here, yielding 
a Q of 25522 ±29. 

In the initial step of the test mass geometry survey, the source and detector masses are 
positioned over the electrostatic shield. The source is tilted on its horizontal axes to ensure 
that the front left corner is the lowest point on the bottom surface. This corner is then used 
as a probe to make a map of the top surface of the detector mass, using a square grid of 
points with a point density of 1.0 mm -2 . The tilt stage of the detector mass is then adjusted 
to optimize the level of the detector, and the surface is mapped again. A similar procedure 
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is then carried out for the electrostatic shield. Finally, a map of the bottom surface of the 
source mass is obtained by touching this surface off against the sapphire hemisphere probe 
fixed to the end of the shield mount. The level is optimized with the source mass tilt stage, 
and the map is repeated. 

In order to obtain the modeshape and amplitude of the driven source mass, the PZT is 
driven at the resonance of the detector and another grid of points is obtained by touching 
the bottom surface of the source mass against the sapphire probe. These measurements are 
done at atmospheric pressure and the touch-offs are determined acoustically. 

The system is then closed and brought down to diffusion pump pressure. The relative 
horizontal positions of the source and detector masses are determined by touching an edge 
of the source against an edge of the detector, and then the horizontal positions are set to 
the fiducial locations. Next, the source is lowered to bring the test masses into momentary 
contact, and their vertical and horizontal positions are recorded. The detector mass is then 
positioned under the electrostatic shield. First it is brought as far forward under the shield 
as possible (without the detector and shield stacks making physical contact), then centered 
between the sides of the shield mount as determined by touch-off signals. The same net 
horizontal translations required to position the detector are applied to the source, bringing 
the test masses into their fiducial positions. 

Finally, the minimum vertical gap is established. The top surface of the detector is 
brought into momentary contact with the bottom surface of the shield, then backed off by 
20 /im. The bottom surface of the source is brought into momentary contact with the top 
surface of the shield, and the source is backed off by 20 /im (the amplitude at the end of the 
source mass) plus an additional 10 /im for safety. 

D. Calibration 

Conversion of the voltage data to observed force is achieved by direct comparison of the 
mean signal to thermal noise. Comparison of the root mean square thermal oscillations of 
the detector mass to the mean displacement of the detector when driven by the source mass 
leads to the expression 
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where V D is the mean of the distribution of voltages in the data sample, Y | V^ 2 " | is the 
component of the standard deviation of the same data sample due to detector thermal 
noise, ks is Boltzmann's constant, and Q,ujq,T, and p d are the detector mechanical quality 



factor, resonant frequency, temperature, and density. Both V D and \j\V T] \ 2 are evaluated 
at a phase of 189°, corresponding to a purely attractive force. The expression ^(f) ■ fir') 
is the projection of the driving force density f(r > ) times the free modal displacement z (r 1 ) 

2 

of the detector at arbitrary point r* on the detector, and z F {f) the mean square free 
modal displacement; both are integrated over the detector volume. A derivation of Eq. |2]is 
given in Appendix A. This expression can be used to convert V D directly to observed force. 
Alternately, the instrument may be calibrated using a calculable electrostatic force or by 
using the reciprocity of the transducer. Both methods agree with calibration based on the 
thermal motion, but the uncertainties involved are substantially larger. 



IV. ANALYSIS 



Generally, the constraints on the Yukawa parameter a for a given range A are obtained 
by ascertaining how large of a Yukawa interaction between the test masses could be present 
and still be consistent with the data in Fig. We present first a simplified analysis based 
on an idealized geometry and neglecting all systematic effects. As the data are consistent 
with a null result limited by detector thermal noise, this preliminary estimate can made by 
comparing the hypothetical Yukawa force between the source and detector with the detector 
thermal noise. 

r-i 

Following the previous analysis of this experiment in Ref. p, the detector modeshape is 
taken to be a pure rotation and the source mass motion a pure translation normal to its 
surface. The source mass is driven at the detector resonant frequency, and only the Fourier 
amplitude of the Yukawa torque at this frequency is effective in driving the detector. From 
Ref. the torque amplitude is 

|AVK)| = 2iTaGp d p s A d RX 2 I l (dz s /X) exp (-g sd /X)[l - exp (-t d /X)][l - exp (-t s /X)\. (3) 

Here, p s is the source mass density, A d is the area of the detector under the source mass 
(half the area of the small rectangle), R is the distance from the edge of the detector to the 
torsion axis (equal to w d /2; see Table HP below). dz s is the source mass amplitude, g sd is the 
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average source-detector gap, I\{dz a /X) is the modified Bessel function, and td and t s are the 
detector and source thicknesses. This expression neglects edge effects and does not include 
the true cantilever modeshape of the source mass. 

The thermal noise torque in the experimental bandwidth on the small detector rectangle 
is found from the mechanical Nyquist expression: 



N 7 



AksT ( mR 2 u \ 



(4) 



\ r \ 3Q )' 

Here, m is the mass of the small rectangle, r is the inverse of the bandwidth, and Q is the 
detector quality factor. 

The signal-to-noise ratio is the ratio of Eq. |3] to Eq. 0] Setting this equal to unity 
and solving for a yields an approximation to the experimental limit curve shown in Fig. 
Evaluating this expression using the mean values of the parameters listed in Tables H] and [H] 
below, assuming an integration time of 77760 s and an average source amplitude equal to 
1/2 the measured tip amplitude, yields a limit in the (a, A) space stronger than the result in 
Fig. El by a small numerical factor for ranges A > 50 /im. The discrepancy increases to about 
two orders of magnitude at A = 5 /im due to the high sensitivity of the Yukawa exponential 
to systematic errors in the source amplitude and the average gap. 

To take the precise geometry and the systematic error into account, the Yukawa force 
between the test masses is computed numerically and constraints on the Yukawa strength 
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a in Fig. [7| are calculated using a maximum likelihood technique [22] . In this approach, the 
interval [a\ , a up ] which contains the true value of a with probability CL is given by 

CL = / p(a\x)da, (5) 

where p(a\x) is the posterior probability density function (p.d.f.) for a given the experi- 
mental data x. The posterior p.d.f. is in turn calculated from Bayes' theorem, 

/ i x L(x\a)it(aL) 

P\ a \ X ) = roo tI I A I A 1 , 1 ( 6 ) 

where L(x\a) is the likelihood function and n(a) is the prior p.d.f. for a. 

In order to account for the effects of various systematics such as test mass geometry, 
density, and other mechanical properties, the likelihood function in Eq. El is replaced by the 
expression 

L'{x\a)= J L(x\a,u)7i(iy)dv, (7) 
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where v represents the set of systematic variables and tt(u) is their prior p.d.f.. In Eq. [?l 
L(x\a, v) is given by 
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where x is the average of the voltages x in the data session, a is their standard deviation, and 
N is the effective number of uncorrelated samples. The term fi(a,v) is the predicted mean 
voltage for a given a and set of systematics if. It is equivalent to V D in Eq. |2l where f(r f ) 
in that equation is understood to be the theoretical value of the Yukawa force density for the 
experimental test mass geometry. A more detailed expression for the likelihood function is 
derived in Appendix B. 



A. Evaluation of the Likelihood Function 

Eq. [7| is evaluated by Monte Carlo integration over the systematics u. For a fixed value 
of the Yukawa range A, 400 points are thrown in the sample space defined by the volume 
dv; the integrand is calculated for each point and added to a running total. 

The parameter x in Eq. |H] represents the mean of the on-resonance data with respect 
to the mean of the off-resonance data (which measures the effective zero), at the phase 
corresponding to the predicted signal. To find x, all data are first projected onto the phase 
of the purely attractive signal. The mean of the off-resonance data is then subtracted from 
that of the on-resonance data. The value of x obtained is —0.44 //V. 

The quantity a/y/N, equal to the standard deviation of the mean of the data, is calculated 
in the following way: First, the projected data sets are partitioned into A part sets of equal 
samples, with A part large enough so that there is very little correlation between partitions. 
The data in each of the partitions are averaged separately and the standard deviation of the 
resulting averages cr av is computed. The standard deviation of the mean is then calculated 
from cr av / J A part . The effective number of uncorrelated samples is the value of N required 
to make a/y/N = a av / ^jN part . Finally, values of a/y/N obtained separately for the on- and 
off-resonance data sets are added in quadrature. The resulting value of a/y/N is 0.82 fiV. 

A list the parameters u needed to evaluate Eq. El is given in Table HI Table |H] lists 
additional systematics which specify the test mass geometry used in the computation of the 
integral / d?r'z F (r') ■ f(r t ) in Eq. [21 
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All parameters in Tables U and |H] are listed with their mean value and error. As the 
information about most of the v is limited, a uniform p.d.f. centered about the mean with 
a width equal to twice the error is assumed for each of the priors tt(u). 

The quantities G, fee, pa, and p s are taken from tabulated values 



J. To obtain the 

thermal noise \J\V T \ 2 , the standard deviation of the off-resonance data (that is, the noise 
due to the JFET amplifier) is subtracted in quadrature from the standard deviation of the 
on-resonance data. The error on \/\V T \ 2 is estimated in the same way as the standard 



deviation of the mean a/yN of the projected data sets. 

Several measurements of the detector mechanical quality factor Q were made immediately 
before and after the data sessions, with the test masses in their fiducial positions. The values 
in Table represent the resulting average and standard deviation of these measurements. 



The quantity y / d 3 f* \z F (f)\ 2 in Eq. |21 the root mean square free modeshape of the 
detector integrated over its volume, is computed numerically from a complete finite element 
model of the detector mass and is treated as a single systematic with a uniform error of 
10%. The free modeshape is normalized to a maximum amplitude of 1 pm. 

For a particular set of the u in Table ITT] a value of the expression / d 3 r'z F (r') ■ f(r r ) in 
Eq. |21 is obtained in a separate Monte Carlo integration. In this calculation, the detector 
curvature, source mass curvature and source mass modeshape are described by second order 
polynomials obtained from fits to the survey data. The detector modeshape is described by 



TABLE I: Systematics in Eq. |2]for evaluation of likelihood function 



Parameter 


Mean 


Error 


Units 


Gravitational constant, G 


6.673 x 1CT 11 


1.0 x 10~ 13 


m kg s 


Boltzmann constant, ks 


1.3806503 x 10~ 23 


2.4 x 10~ 29 


JK- 1 


Detector density (tungsten), pd 


1.93 x 10 4 


1.9 x 10 3 


kg m -3 


Source density (tungsten), p s 


1.93 x 10 4 


1.9 x 10 3 


kg m -3 


Thermal noise voltage, \J\V T \ 2 


6.09 x 10~ 5 


2.3 x 10~ 6 


V 


Mechanical quality factor, Q 


2.5522 x 10 4 


29 


(NA) 


Resonant frequency, ujo/2it 


1173.085 


0.015 


Hz 


Temperature, T 


305.0 


0.1 


K 


Integrated rms free modeshape, 


sjfd^l z F (f)\ 2 5.87 xlO- 11 


5.9 x 10- 12 


m 5/2 
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TABLE II: Systematics for evaluation of / d 3 r'z F (r') ■ f(r / ). All units are meters. See text for 
definition of parameters. 

Parameter Mean Error 

Detector length, l d 5.0800 x 10~ 3 6.4 x 10~ 6 

Detector width, w d 1.14550 x 10~ 2 6.4 x 10~ 6 

Detector thickness, t d 1.950 x 10" 4 6.4 x 10~ 6 

Source width, w s 7.0000 x 10~ 3 6.4 x 10~ 6 

Source thickness, t s 3.048 x 10~ 4 6.4 x 10~ 6 

Touch gap, g sd 1.080 x 10~ 4 6.4 x 10~ 6 

Source amplitude, dz s 1.87 x 10~ 5 3.2 x 10~ 6 



a similar function obtained from the finite element model. The curvature and modeshape 
functions and the v are used to calculate bounding boxes for the source mass and the small 
rectangle of the detector mass. For a particular position of the source mass, 1 x 10 5 point 
pairs are thrown. For those pairs landing within the bounding boxes a force density f [f") is 
calculated and added to a running total. A total of 30 different phase positions of the source 
mass between and 2n radians are sampled, with the range defined by the amplitude dz s . 

The vertical separation of the test masses is computed from the touch gap, the test 
mass thicknesses (t d ,t s ), and a correction calculated from the functions describing the test 
mass surface curvature. The touch gap (g sd ) is the vertical difference in test mass positions 
between where the opposing test mass surfaces touch and where they are in their fiducial 
positions. The error on these quantities is taken to be twice the resolution of the translation 
stages, or 6.4 /mi, as two measurements are needed to determine each. This is also the 
case for the errors in the test mass widths (w d ,w s ) and detector length (l d ). An exception 
is the source mass tip amplitude (dz s ), which is the average of five measurements. Also, 
the accuracy with which the tip of the source mass is aligned with the back edge of the 
small detector rectangle is 320 microns, as it is aligned by eye and taken to have an error 
corresponding to a complete turn of the translation stage screw. The accuracy with which 
the long edge of the source mass is aligned to the torsion axis is the usual 6.4 microns, as 
this position is determined from touch-offs. All geometry measurements are repeatable over 
several weeks to within the resolution of the translation stages. 
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FIG. 6: Likelihood function (Eq. [7J) dependence on a for the case A = 20 /im. 

Fig. |U1 shows the likelihood function (Eq. [7J) as a function of a for the case A = 20 /im, 
a range at which the experiment has good sensitivity relative to previous experiments. The 
factor (o^/2ixej has been suppressed as it has no bearing on the posterior p.d.f. for 
a (Eq. IHJ). The maximum is shifted slightly toward negative (repulsive) values of a, in 
accordance with the offset of the mean of the on-resonance data with respect to the off- 
resonance data (Fig. EJ). 

B. Limits on Yukawa Parameters 

The function in Fig. El is substituted into Eq. |U1 and integrated numerically (Eq. |3J) over 
an interval [a\ , a up ] to obtain the limit on a. The interval is centered around the maximum 
of the likelihood function and the length adjusted by trial and error to obtain a confidence 
level of 95%. In Eq. El a uniform prior distribution for a is used. For the case A = 20 /im, 
the resulting 95% confidence interval is [—5.60 x 10 3 < a < 3.56 x 10 3 ]. 

The likelihood function is recomputed and the integration repeated for several values 
of A between 5 /im and 500 /im. The resulting 95% confidence level limit curve in the 
(a, A) parameter space is shown in Fig. [7| For each value of A, the slightly weaker limit 
corresponding to a repulsive interaction is shown. The new limit is close to three orders 
of magnitude more sensitive than the previous experimental limits in the range between 10 
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and 100 /zm. 
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FIG. 7: 95% confidence level limit on the Yukawa strength a as a function of the range A for 
this experiment (bold solid curve), together with limits from previous experiments (bold dashed 
curves) and theoretical predictions (fine lines). 



V. CONCLUSION 

An experimental search for new macroscopic forces at short distances has been conducted 
using 1 kHz mechanical oscillators as test masses with a stiff conducting shield between them. 
No evidence for a resonant signal above detector thermal noise has been observed. Based 
on these measurements, roughly three new square decades in the (a, A) parameter space are 
ruled out in the range 10 /im < A < 100 /im. About half of the remaining parameter space 
for the gluon and strange modulus forces is now excluded. Our limit crosses the predicted 
line for the dilaton at 36 microns, corresponding to a lower limit on the dilaton mass of 
5.5 x 10~ 3 eV. 
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From Fig. [71 the sensitivity of the experiment maximizes at about 4 times gravitational 
strength above A = 200 /im. For the planar test mass geometry of the current experiment, 
the Newtonian signal is roughly an order of magnitude below the thermal noise at this 
range. Detection of this signal should be possible with more statistics and an optimized 
geometry. The small test masses in this experiment allow for relatively easy control of the 
minimum separation, currently limited to 100 fim by the thickness of the electrostatic shield. 
Preliminary tests with stretched metal membranes suggest that shields of sufficiently low 
compliance can be made with thicknesses as small as 10 /im. Work is underway toward 
an experiment with a stretched membrane shield in place of the sapphire plate. If the 
backgrounds can be controlled, this experiment could improve limits between 10 and 50 /im 
by at least another order of magnitude. 
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APPENDIX A: CALIBRATION 

The mean signal from the lock-in amplifier data can be converted to observed force by 
direct comparison to the detector thermal noise. The mean square thermal displacement of 
the detector mass ^(r)) 2 at the position f of the capacitive transducer probe is related to 
the temperature T via the equipartition theorem: 

l -m{r)WW^ = \k B T. (Al) 

Here, ujq is the resonant frequency of the 5th detector mode and m(r) is the modal mass of 
the detector oscillator at the position of the probe. 

For a distributed mass with free modal displacement z F (r') at an arbitrary point r* driven 
on resonance, the driven displacement amplitude dz ^) due to a differential force amplitude 



21 



dF(r l ) is given by 



dz D (7) = -j 



Q 



■z F (f* 



z (f 1 ) ■ dF(f) 



(A2) 



where Q is the detector mechanical quality factor and the dot product picks out the compo- 
nent of the force parallel to the modal displacement. (Note the use of e +JaJ * time dependence. 
The differential amplitude at the probe is 

z F (r) 



dz D {r) 



dz D (?\ 



(A3) 



Substituting Eq. IA3I into Eq. IA2I and integrating over r* yields the amplitude of the detector 
at the position of the transducer due the total force on the detector: 



(A4) 



Jd 3 P\z F (P)\ 2 

Here the force amplitude density f(r ! ) is defined by f{f) = dF(r , )/d 3 f f , pa is the detector 
mass density, and the definition of the modal mass 



m(f) 

has been used. 

Combining Eqs. IA11 IA41 and IA5I yields 

£ D (r) Q 



pdl 


Z F (F) 






z F {f)\ 2 



(A5) 



z F {r) fd 3 7z F {f*) ■ /(f) 



^(f)! 2 juj y/k B Tp d \z F {r) 



fd?r > \z F {r > )\' 



(A6) 



The linear response of the capacitive transducer insures that the ratio z D (r)/\/\ : ' (r) 
related to the measured voltages on the lock-in amplifier by 



is 



z D (r) 



V D 



^(rir 



\V T \' 



(A7) 



where V D is the mean of the distribution of voltages in the data sample and Vl^ T | 2 * s the 
component of the standard deviation of the same data sample due to detector thermal noise. 
Both quantities are evaluated at the phase of interest. Substitution of Eq. IA7I into Eq. IA6I 
yields Eq. H 
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APPENDIX B: LIKELIHOOD FUNCTION 



In addition to the unknown parameter a corresponding to the Yukawa strength, the 
signal depends on a set of systematics v which includes test mass geometry, density, detector 
temperature and mechanical properties. The effects of uncertainties in these quantities can 
be accounted for in the posterior p.d.f. (Eq. (BJ) by integrating over the systematics: 



p(a\x) = J p(a,u\x)du. 



(Bl) 



Assuming the prior joint p.d.f. for a and u factorizes, this is equivalent to replacing the 
likelihood function L(x\a) with 



L'{x\a) = J L(x\a,is)n(is)dis. 



(B2) 



Here, 7t(v) is the prior p.d.f. of the v. The function L(x\a,v) is the joint p.d.f. for the 
data, regarded as a function of both the unknown Yukawa strength a and the systematics 
v. 

L(x\a,u) = Y[F(xi\a,v). (B3) 

i 

For the present experiment, the data X{ are the sampled voltages from the lock-in am- 
plifier, which are Gaussian distributed (Fig. |1J). Therefore the function F is taken to be 
the Gaussian distribution G(xi\/j,(a,u),cr), where fi(a,u) is the predicted mean voltage for 
a given a and set of systematics u, and a is the observed standard deviation: 

( 1 \ 



F(xi\a, v) 



exp 



X; 



(B4) 



The product is 



N 



Y[F(xi\a,v) 



' 1 



exp 



N 



N 



exp 



(x - n(a, v)f 



2(<t/VN) 



(B5) 



where x is the average of the voltages Xj and N is the effective number of uncorrelated 
samples (the value of N such that a /y~N is the standard deviation of the mean of the actual 
data). 

The predicted mean fi(a,u) can be replaced by Eq. |2]for V D if /(f) is identified with 
the theoretical value of the Yukawa force density for a test mass geometry corresponding to 
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a particular set of parameters v. As a Yukawa force density, f(r r ) is linear in the parameter 
a, the gravitational constant G, and the source mass and detector mass densities p s and pd- 
f(r > ) = aGpspdif). Substituting into Eq. Eland factoring these terms out of the integral 
yields 

v \uj Vk B TJ . s Jfd3r'\z F (r>)\ 2 

Substituting Eq. IB6I into Eq. IB5I and the result into Eq. IB2I yields the likelihood function 
L'(x\a). 
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